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Abstract. The production of K + mesons in proton-nucleus collisions from 1.0 to 2.3 GeV is analyzed 
with respect to one-step nucleon-nucleon (NN — > NYK + ) and two-step zi- nucleon (AN — > K + YN) 
or pion-nucleon (nN — > K + Y) production channels on the basis of a coupled-channel transport approach 
(CBUU) including the kaon final-state-interactions (FSI). Momentum-dependent potentials for the nucleon, 
hyperon and kaon in the final state are included as well as K + elastic rescattering in the target nucleus. 
The transport calculations are compared to the experimental K + spectra taken at COSY-Jiilich. Our 
systematic analysis of K + spectra from 12 C, 63 Cu, 107 Ag and 197 Au targets as well as their momentum 
differential ratios gives a repulsive K + potential of 20 ± 5 MeV at normal nuclear matter density. 



PACS. 13.60.Le Meson production - 13.75.Jz Kaon-baryon interactions 
- 24.40.-h Nucleon-induced reactions 



14.40.Aq Pi, K, and eta mesons 



1 Introduction 

The production of heavy mesons in p+A reactions at bom- 
barding energies far below and close to the free nucleon- 
nucleon threshold is of specific interest PQ-^ZI as one hopes 
to learn either about cooperative nuclear phenomena and/ 
1 or about high momentum components of the nuclear many- 
body wave function that arise from nucleon-nucleon cor- 
relations. Especially K + mesons have been considered as 
promising hadronic probes [T%lll9ll20| due to the rather 
moderate final state interaction, which is a consequence 
of strangeness conservation. 

Experiments on production from nucleus-nucleus 
collisions at SIS energies of 1- 2 A- GeV have shown that 
in-medium properties of the kaons are seen in the collec- 
tive flow pattern of K + mesons both, in-plane and out-of- 
plane, as well as in the abundancy of antikaons [2"Tll22ll23) . 
Thus in-medium modifications of the mesons have become 
a topic of substantial interest in the last decade triggered 
in part by the early suggestion of Brown and Rho [21] , that 
the modifications of hadron masses should scale with the 
scalar quark condensate < qq > at finite baryon density. 

As demonstrated in the pioneering work of Kaplan and 
Nelson j23l2*5| kaons and antikaons couple attractively to 
the scalar nucleon density with a strength proportional 
to the KN — £ constant Skn , which is not well known 
at present and may vary from 270 to 450 MeV. Further- 
more, a vector coupling to the quark 4-current - for van- 
ishing spatial components - leads to a repulsive potential 
term for the kaons and to an additional attractive term for 
the antikaons. In fact, the recent observation of a rather 



narrow strange tribaryon state with mass 3.115 GeV in 
Ref. 23 ~ observed in the 4 He(K~ , p)X reaction with an- 
tikaons at rest - has given further experimental evidence 
for strong attractive K~ potentials in (light) nuclei. 

Nevertheless, the actual kaon and antikaon self ener- 
gies (or potentials) at normal nuclear matter density (and 
above) are quite a matter of debate - due to higher or- 
der terms in the chiral expansion - especially for the an- 
tikaon 28 29, 30' • Moreover the momentum-dependence of 
their self energies is widely unknown (except for a disper- 
sion analysis in Ref. |31| ) since most Lagrangian mod- 
els are restricted to s-wave interactions or only include 
additional p-waves. It is thus mandatory to perform ex- 
perimental studies of the (anti-) kaon properties under 
well controlled conditions, e.g. in proton-nucleus reactions, 
where one probes the (anti-) kaon self energies at normal 
nuclear matter density po sa 0.16 fm -3 and below. Fur- 
thermore, by gating on kaon momenta in the laboratory, 
one might be able to obtain information on the momentum 
dependence of the self energies, too. 

K + production in p+A collisions at subthreshold ener- 
gies has been observed experimentally more than a decade 
ago by Koptev et al. ^H] at bombarding energies from 0.8 
to 1.0 GeV. However, only total K + yields could be ex- 
tracted at that time at subthreshold energies. In more 
recent years differential K + spectra have been measured 
down to 1.2 GeV for 12 C targets at 40 [32] (SATURNE) 
or 90 in the laboratory [33] (CELSIUS). Unfortunately, 
the different experiments have no overlap in acceptance 
and the interpretation of the data, if compatible at all, 
remains vague (14|34j. First data on the full momentum 
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distribution at forward angles have been presented by the 
ANKE Collaboration at COSY-Jiilich gHl for K+ mesons 
from p + 12 C reactions at 1.0 GcV in 2001 ^ and at 
higher bombarding energies in the last year |37| . 

In this study we use the coupled-channel (CBUU) trans- 
port model that has been developed first in Ref. J3B] for 
the description of nucleus-nucleus collisions and later on 
employed for the simulation of pion- and proton-nucleus 
reactions [3*Tll39ll40| . too. For applications to K ± produc- 
tion in nucleus-nucleus collisions at SIS energies we re- 
fer the reader to Refs. In this model the effects of 
momentum-dependent self energies for all hadrons can be 
studied explicitly as well as their production and propa- 
gation in the nuclear medium. The actual version of the 
CBUU transport model employed in our present analysis 
is identical to that described in Ref. . We thus discard 
an explicit description and refer the reader to Ref. for 
technical details and the explicit momentum-dependent 
potentials employed for the nucleon and hyperon degrees 
of freedom. 

The paper is organized as follows: We very briefly re- 
capitulate the concepts of the CBUU model in Section 
2 and step on with the actual results for K + spectra in 
p + A collisions at COSY energies in Section 3 in com- 
parison to the data of the ANKE collaboration In 
Section 4 we compare the calculated ratios of K + spectra 
for different kaon potentials with the corresponding data 
at various bombarding energies |37| . This systematic ana- 
lysis will lead to an approximate extraction of the strength 
of the K + potential in the nuclear medium. A summary 
concludes this paper in Section 5. 



2 Brief reminder of the CBUU transport 
model 

In this work we perform our analysis along the line of the 
CBUU 1 approach [3H1 which is based on a coupled set 
of transport equations for the phase-space distributions 
fh(x,p) of hadron h, i.e. fJ21 

(7T M - n v dfJJ v h ~ MZ8>U%) d£f h (x, p) (1) 
+ (n v d£UZ + M* h d^) d»f h {x,p) 

= J d2d3dA[G^G] 12 ^M... 

x <5 4 (i7 + n 2 -n 3 - n 4 ) 

x {fh a (X, P3)fh 4 (X, Pi)fh (X, p)fh 2 (X, P2) 

~ fh(x,p)f fl2 (x,P2)fh 3 (x,p 3 )fh 4 (x,Pa)} • (2) 

In Eq. Q Ufr(x,p) and U£(x,p) denote the real part of the 
scalar and vector hadron self energies, respectively, while 
[G*G] 12-34... S 4 (n + 77 2 - i7 3 - 7T 4 . . .) is the 'transition 
rate' for the process 1 + 2^3 + 4 + ... which is taken to 
be on-shell in the semiclassical limit adopted. The hadron 

Coupled-Channel Boltzmann-Uehling-Uhlenbeck 



quasi-particle properties in |TjQ) are defined via the mass- 
shell constraint 

5(n,n»-M* h 2 ) , (3) 

with effective masses and momenta (in local Thomas- Fermi 
approximation) given by |42j 

Ml{x,p) = M h + U%(x,p) 
n»{x,p)=p> i -U£{x,p) , (4) 

while the phase-space factors 

fh(x,p) = 1 ± fh(x,p) (5) 

are responsible for fermion Pauli-blocking or Bose enhance- 
ment, respectively, depending on the type of hadron in the 
final/initial channel. The transport approach Q is fully 
specified by U^{x,p) and U£(x,p) (// = 0,1,2,3), which 
determine the mean-field propagation of the hadrons, and 
by the transition rates G^G 5 4 (. . .) in the collision term 
(5), that describes the scattering and hadron production/- 
absorption rates. 

The scalar and vector mean fields U% and U£ for nu- 
cleons are modeled in line with Ref. and presented 
explicitly in Fig. 1 of Ref. 34 for the resulting real part 
of the nucleon potential Un{p)- We recall that this nu- 
cleon potential (at density po) is attractive for momenta p 
below 0.6 GeV/c and repulsive above. Apart from the nu- 
clear potentials each charged hadron additionally moves 
in the background of the Coulomb potential that is gener- 
ated by the charged hadrons themselves. In case of proton- 
nucleus reactions - with the nucleus at rest - this essen- 
tially amounts to a Coulomb acceleration in the final state 
for positively charged hadrons and a deceleration for nega- 
tively charged particles. Note, that for heavy nuclei like Pb 
or Au the Coulomb potential in the nuclear interior (for 
p, tt + ,K + ) is ~ 20 MeV, i.e. of the same order of magni- 
tude as the 'expected' repulsive K + nuclear potential [341 

EE]. 

The hyperon mean fields, furthermore, are assumed 
to be 2/3 of the nucleon potentials. In the present ap- 
proach, apart from nucleons, A, AT(1440), iV(1535), A, £ 
with their isospin degrees of freedom, we propagate ex- 
plicitly pions, K + , and 77's and assume that the pions as 
Goldstone bosons do not change their properties in the 
medium. The kaon potential is taken as momentum inde- 
pendent (cf. Ref. 3 1 i . while its strength is varied in the 
dynamical calculations from to 20 MeV at nuclear mat- 
ter density po assuming a linear increase with the actual 
nucleon density p(r). 

The calculation of 'subthreshold' particle production 
is treated perturbatively in the energy regime of inter- 
est due to the small cross sections involved. Since we 
work within the parallel ensemble algorithm in the CBUU 
approach, each parallel run of the transport calculation 
can be considered approximately as an individual reaction 
event, where binary reactions in the entrance channel at 
given invariant energy y/s lead to final states with 2 (e.g. 
K + Y in irB channels), 3 (e.g. for K + YN channels in BB 
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collisions) or 4 particles (e.g. KKNN in BB collisions) 
with a relative weight Pj for each event i which is defined 
by the ratio of the production cross section to the total 
hadron-hadron cross section. We thus dynamically gate on 
all events where a K + meson is produced initially. Each 
strange hadron production event i is represented by ~ 10 3 
testparticles for the final strange hadron j with individual 
weight Wj such that the sum of the weights Wj over j re- 
produces the individual production probability Pi and the 
distribution in momenta (multiplied with the NN or itN 
cross section) describes the differential production cross 
section (cf. Ref. [S] for details). 

For the present study the production of pions by pN 
collisions in p + A reactions as well as the total kaon cross 
sections in pN and ir N collisions are of relevance. The pion 
production cross section from NN interactions is based on 
the paramctrization of the experimental data by Ver West 
and Arndt and implemented in the transport model 
as described in Ref. The cross sections for the chan- 
nels irN — > KY, where Y stands for a hyperon A, E, are 
taken from the analysis of Huang et al. @B] and essentially 
correspond to the experimental data for the different itN 
channels in vacuum (or 'free' space). Note, that in addi- 
tion to the early studies in |16U17| the channels with a 
£ hyperon are taken into account. All cross sections are 
reparametrized as a function of the invariant energy above 
threshold \/s — y/so where ^/so denotes the threshold 
for the individual channel given by the sum of the hadron 
masses in the final state of the reaction. We recall that 
the differential cross sections for binary reactions are fully 
determined by the total cross sections when assuming S- 
wave dominance. 

The question arises, furthermore, about the in-medium 
production cross sections - essentially at density p n " if P°~ 
tentials or self energies are involved. Here we employ the 
assumption that the production matrix element squared 
\M\ 2 does not change in the medium and the change of 
cross section can be described by a change of the available 
phase space. This notion is guided by the experimental 
observation that meson production from pp collisions is 
well described by phase space 07J if final state interac- 
tions (FSI) between the hadrons are taken into account. 
These FSI, furthermore, are found to be dominated by 
the final baryon-baryon interaction, which is very strong 
in the pp,pn or pA channels in free space. On the other 
hand, such FSI's are screened in the nuclear medium due 
to the surrounding nucleons such that in-medium produc- 
tion cross sections for mesons are expected to vary essen- 
tially with the available phase space. It is thus sufficient 
to shift the threshold in a pN collision to 

y/4 = n N (p N ) + n A ( PA ) + n K (p K ) (6) 

using Q), where the momenta pj denote the relative mo- 
menta with respect to the nuclear matter rest frame and 
IIj(p) is the energy of hadron j including the potentials 

(4). ' 



3 Comparison to experimental spectra 

In this Section we show the influence of nucleon and kaon 
potentials on K + spectra from p + 12 C and p + 197 Au re- 
actions and compare our calculations to the experimental 
K + spectra available from 1.0 GcV — 2.3 GeV bombard- 
ing energy from the ANKE collaboration 37 that has 
taken K + spectra in forward direction for 9i a b < 12° up 
to momenta of ~ 500 MeV/c. 




P K + (MeV/c) P K + (MeV/c) 

Fig. 1. The calculated differential K + spectra for p + 12 C 
from 1.0 to 2.3 GeV for 6i ab < 12°. The dotted line is ob- 
tained from CBUU calculations without baryon and kaon po- 
tentials, the dashed line shows the results with baryon poten- 
tials included while the solid line corresponds to calculations 
including both, nucleon and kaon potentials (+20 MeV at nu- 
clear density po). The shaded areas indicates the contributions 
from the two-step mechanisms pN — > AN, AN — > K + YN and 
pN — > tvNN, ttN — » K + Y , respectively, for the case of nucleon 
and kaon potentials, such that the difference to the solid line 
corresponds to the primary pN production channel. 

The calculated differential K + spectra for p+ 12 C from 
1.0 GeV to 2.3 GeV for 6i ab < 12° are displayed in Fig. 
n and for p + 197 Au in Fig. [21 The dotted lines are ob- 
tained from CBUU calculations without baryon and kaon 
potentials, the dashed lines show the results with baryon 
potentials included while the solid lines correspond to cal- 
culations with both, nucleon and kaon potentials. In the 
latter case we have assumed a kaon potential of +20 MeV 
at density po. At the low bombarding energy of 1.0 GeV 
the nuclear potential in the final state is dominantly at- 
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Fig. 2. The same as Fig. 1 for p + Au reactions from 1.0 to 
2.3 GeV bombarding energy. 



tractive and we obtain an enhancement of the K + yield by 
about a factor of 2 - relative to the calculations without 
nuclear potential - whereas the additional repulsive K + 
potential leads again to a decrease by a factor ~ 2-3. For 
1.2 GeV the net effect of the nucleon potential is close to 
zero; its repulsion in the final state for higher bombarding 
energies becomes visible in a reduction of the K + spec- 
trum relative to the 'free' case. On the other hand the 
effect of the repulsive kaon potential is always a reduc- 
tion of the K + spectra, which is most pronounced at very 
low bombarding energies - when comparing the solid lines 
relative to the dashed lines - and becomes small at high 
bombarding energies. 

The shaded areas in Figs. ^ an d 121 indicate the contri- 
butions from the two-step mechanisms pN — > AN, AN — ► 
K+YN and pN -> ttNN,ttN -> K+Y, respectively, for 
the case of nucleon and kaon potentials included (solid 
line). Thus the role of secondary (A and pion induced) 
reaction channels is clearly visible from Figs. 1 and 2 by 
comparing the shaded area to the solid line, that corre- 
spond to the total spectra. At T/ Q b = 1.0 GeV the domi- 
nant fraction of the K + yield is due to the secondary chan- 
nels in line with the earlier calculations in Refs. |16II17| . 
Consequently, one does not probe high momentum com- 
ponents of the nuclear wave function by K + spectra in a 
direct way. However, with increasing bombarding energy 
the role of the two-step reactions diminishes gradually. At 
Tiab = 2.3 GeV the secondary channels in case of a Au 
target amount to about 30% and in case of a C target 
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Fig. 3. Comparison of the K + spectra from the CBUU cal- 
culations (solid lines) with the data from the ANKE collabo- 
ration |3fcil37| for p + 12 C. The calculations include both the 
momentum-dependent nucleon potential as well as a repulsive 
K + potential of +20 MeV at density po- 



to less than 20%. This relative change with target mass 
number is attributed to the fact that for the small 12 C 
target only a fraction of the high energy pions rescatters 
in the target and produces K + Y pairs. Moreover, the role 
of the secondary channels decreases with increasing kaon 
momentum such that the high momentum tail of the K + 
spectra is dominated by the first chance pN production 
channel in line with the studies by Paryev pEj . 

How do these calculations compare to the data from 
the ANKE collaboration 36,37 and what can we learn 
from that? To answer this question we show in Fig.|21our 
calculations that include the baryon and K + potentials 
(solid lines) with the measured spectra for p + 12 C from 
Refs. |36U37j . Whereas the description of the experimen- 
tal spectrum looks almost perfect at 1.0 GeV the data 
are slightly underestimated at Ti a t — 1.2 and 1.5 GeV 
in the maximum. At these energies the maximum of the 
spectrum is better described without kaon and nuclear 
potentials, but then the spectra no longer match in the 
low momentum part. Thus a repulsive K + potential is 
needed to properly describe the shape of the spectrum es- 
pecially for low momenta. On the other hand, channels 
such as pN% — > d + tt followed by a secondary reaction 
7r + N2 — > K + + Y are not included in the present trans- 
port calculations. Such channels have been found in Ref. 
|49j to contribute to the kaon spectrum at bombarding 
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energies ~ 1.2 GeV but to be insignificant at higher en- 
ergies of 2.0 to 2.3 GeV. As can be seen from Fig. the 
calculations - including nucleon and kaon potentials - are 
in a good agreement with the experimental spectra again 
at 2.0 and 2.3 GeV which is in line with the findings of 
Ref. HI]. 
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Fig. 4. The ratio of the calculated differential K + spectra 
from p + Au to p + C reactions at 1.5, 1.75 and 2.3 GeV in 
comparison to the data from Ref. |37|. The different solid lines 
correspond to calculations employing different kaon potentials 
at normal nuclear density po- Starting from the left the results 
correspond to repulsive kaon potentials of 0, 10 and 20 MeV 
at po- 



4 K + spectral ratios 

The differential K + momentum spectra in Figs. 1 to 3 
differ in their low momentum part - at fixed bombarding 
energy - when varying the target mass number. In Refs. 
|34U44j it has been pointed out that the variation in the 
low K + momentum spectrum is a consequence of i) the 
different acceleration of the charged kaon in the Coulomb 
field of the target and ii) an additional acceleration in the 
repulsive kaon potential in the nuclear medium (target 
nucleus). The combined effect of both the Coulomb and 
nuclear repulsion is most easily seen when comparing the 
K + spectra from different targets. To this aim we show in 
Fig.^Jthe ratio of the spectra from Au to C targets at 1.5, 
1.75 and 2.3 GeV, where this effect is most pronounced (cf. 
Ref. 01]). The different solid lines in Fig. ^correspond to 
calculations employing different nuclear kaon potentials of 
magnitude 0, 10 and 20 MeV at normal nuclear density 
po (starting from the left). 

It is seen that the ratio for (Au/C) is approximately 
independent on the bombarding energy and shows a max- 
imum at about 230 MeV/c. When neglecting kaon repul- 
sion from nuclear forces - including only Coulomb repul- 
sion - this maximum shows up slightly below 200 MeV/c, 
which is clearly too low in comparison to the data (cf. |37l 
I44|h We recall that the relative strength of the momen- 
tum shift in the forward K + spectra is proportional to the 
square root of the sum of Coulomb and nuclear potentials, 

Ap a y/2M K {Ucoui + U K ), (7) 

which implies that an additional repulsive nuclear K + po- 
tential is needed to describe the data. Since the strength 
of this potential - linear in the nuclear density at the K + 
production point - is a 'parameter' in the transport calcu- 
lations, we may in turn determine the strength of this po- 
tential in comparison to the data. A x 2_ fit gives Uk(po) ~ 
20 ± 5 MeV for the ratio of p + Au to p + C reactions. 

We mention that the transport calculations have a 
limited accuracy due to low statistics in some regions of 
phase space. This holds true especially for the low mo- 
mentum spectrum of the kaons in the acceptance of the 
ANKE spectrometer and has an impact also on the ra- 
tio of the spectra from different targets. In the x 2_nt ad- 
dressed above the statistical inaccuracy of the calculations 
has been included. 

It is of further interest whether the interpretation of 
the momentum differential kaon spectrum - discussed above 
- is compatible with the ratios from other targets, too. To 
this aim we show in Fig. the same ratios as in Fig. 0] 
but for Ag/C and Cu/C at 2.3 GeV in comparison to the 
data from Ref. [37] . It is clearly seen that the experimen- 
tal ratios are consistent with the transport calculations for 
Uk(po) = +20 MeV though the sensitivity to the strength 
of the kaon potential is reduced. This is basically due to 
the fact that the surface area - relative to the total vol- 
ume - of the lighter targets is larger than for the heavy 
Au nucleus. 

Nevertheless, the agreement between the CBUU calcu- 
lations and the data is remarkable and demonstrates that 
the underlying K + dynamics is rather well understood. 
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Fig. 5. The same ratios as in Fig. 0]but for Ag/C and Cu/C 
at 2.3 GeV in comparison to the data from Ref. |37) . 

5 Summary 

In this work we have studied the production of K + mesons 
in proton-nucleus collisions from 1.0 to 2.3 GeV with re- 
spect to one-step nucleon-nucleon and two-step Z\-nucleon 
or pion-nucleon production channels on the basis of a 
coupled-channel transport approach (CBUU) including dif- 
ferential transition probabilities from 7riV reactions (cf. 
Ref. [S3)- We have included the kaon final state interac- 
tions, which are important for heavy targets, and incor- 
porated the effects of momentum-dependent potentials for 
the nucleon, hyperon and kaon in the nucleus. A compar- 
ison of the calculations to the experimental K + spectra 
taken at COSY-Jiilich has shown that the momentum dif- 
ferential K + yield can be reproduced in magnitude and 
shape within an accuracy better than 30% at all bom- 
barding energies. 

The detailed calculations have demonstrated that the 
spectra show a substantial sensitivity to the potentials 
and their momentum dependence. At low bombarding en- 
ergies of ~ 1.0 GeV the net attractive potentials for the 
nucleon and the /1-hyperon in the final state lead to a 
relative enhancement of the K + spectra while at higher 
bombarding energies (~ 2 GeV) the baryon potentials are 
repulsive and thus suppress K + production relative to the 
free case. This phenomenon is seen in the excitation func- 
tion of the K + cross section when varying Ti a b from 1.0 - 
2.3 GeV (cf. Fig. 3). Furthermore, the shape of the spec- 
trum for low K + momenta in the laboratory is very sensi- 
tive to both, Coulomb and nuclear kaon potentials, since 



the kaons are accelerated by both forces when leaving the 
nuclear environment and propagating to the continuum. 
The relative strength of this momentum shift in the for- 
ward K + spectra is proportional to the square root of the 
sum of both potentials 0. Thus the K + spectral shape 
at low momenta (or kinetic energies Tk) allows to deter- 
mine the strength of the K + potential from experimental 
data in an almost model independent way especially when 
comparing kaon spectra from light and heavy targets at 
the same bombarding energy (cf. Refs. [3711441150) 1. In fact, 
the shape of such ratios is practically independent of bom- 
barding energy but sensitive to the strength of the kaon 
repulsion in the medium. Our systematic study has given 
a value of 20 ± 5 MeV at normal nuclear density po which 
agrees with the result from the scattering length approxi- 
mation 

We point out that our calculations slightly underpre- 
dict the spectra at 1.2 and 1.5 GeV indicating that ad- 
ditional reaction channels should be at work which have 
not been included in the CBUU calculations. In fact, the 
experimental studies by the ANKE collaboration in Ref. 
0Hj have shown that there is an additional K + production 
from the secondary channel p + N\ — » d + it; ir + N2 — > 
K + + A at a bombarding energy of 1.2 GeV in p + C 
reactions. However, no signals for this channel have been 
seen at higher bombarding energies of ~ 2 GeV. This in- 
terpretation is well in line with our present transport cal- 
culations since the description of the available spectra at 
2.0 and 2.3 GeV for p + 12 C implies that further reac- 
tion channels do not contribute significantly at the higher 
bombarding energy. It remains presently an open question 
whether a cooperative production mechanism by two- or 
three-nucleon clusters - not included in our transport cal- 
culations - might contribute and be identified from the 
experimental side in p + A reactions at very low energies. 

The authors like to acknowledge valuable discussions 
with M. Buschcr, V. Koptev, M. Nekipelov, P. Senger, A. 
Sibirtsev, H. Stroher and C. Wilkin on various issues of 
this study. 
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